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Abstract 

We establish necessary and sufficient conditions for the shock statis- 
tics to approach self-similar form in Burgers turbulence with Levy pro- 
cess initial data. The proof relies upon an elegant closure theorem 
of Bertoin and Carraro and Duchon that reduces the study of shock 
statistics to Smoluchowski's coagulation equation with additive kernel, 
and upon our previous characterization of the domains of attraction of 
self-similar solutions for this equation. 
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1 Introduction 



The construction of stochastic processes that are also weak solutions to 
the equations of fluid mechanics is one approach to rigorous mathematical 
theories of turbulence. This is poorly understood at present, and we must 
settle for insights from vastly simplified model problems. We consider the 
invisicid Burgers equation 

d t u + d x (^- S j=0, t>0,x£R, u(x,0) = u o (x), (1) 

with random initial data no- The problem is to determine the statistical 
properties of the Cole-Hopf (entropy) solution u(x,t) to (1), given the sta- 
tistical properties of uq. There is a large literature on the subject; we refer 
to Burgers' book [10] and the more recent survey articles [17, 25, 36]. The 
problem was proposed by Burgers as a model for turbulence in incompress- 
ible fluids, but it has several well-known flaws in this regard. 

Explicit solutions play a special role in the theory. Burgers studied the 
case when no is white noise in his monograph [10]. His work remains the 
foundation for several rigorous results, which culminate with the complete 
solution by Frachebourg and Martin for the velocity and shock statistics (see 
[19] and references therein). The case when uq is a Brownian motion has 
attracted much attention since the work of She, Aurell and Frisch [34] and 
Sinai [35]. An elegant solution to this problem was obtained by Bertoin [5] 
and Carraro and Duchon [12]. More generally, these authors considered 
initial data that comprise a Levy process with only downward jumps (i.e., 
shocks) . A Levy process X x (x > 0) is a continuous-time random walk with 
stationary and independent increments. It is determined completely by its 
characteristic exponent Vf, satisfying ~E(e lkXx ) = e ~ x ^( k ) } via the celebrated 
Levy-Khintchine formula 

*(A;) =ibk + — - + / (l-e lks + iksl lsl<1 )U(ds), k£R. (2) 

The process X is the superposition of three independent processes related 
to this formula: a Brownian motion with variance a 2 and drift b £ R, a 
compound Poisson process with jump measure IH| s |>i5 an d a pure jump 
martingale with jump measure ni| s | <:L (see [4, Ch.l]). The measure LT is 
arbitrary, subject to the condition A s 2 )H(ds) < oo, where a A b means 
min(a, b). 
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We say X is spectrally negative if II is concentrated on the half-line s < 0. 
In all that follows we assume 



As we show below, we can always reduce to the case where uo(x) has zero 
mean K(uo(x)) for all x, and J° (|s| A s 2 )H(ds) < oo. It is then more 
convenient to use the Laplace exponent 



where A((a,oo)) = n((-oo, -s)) for every s > 0, so that E(e qX *) = e x ^ 
and E(X X ) = x^'(0) = 0. 

For this class of initial data, Bertoin proved a remarkable closure prop- 
erty for the entropy solution of (1), namely: x t— > u(x,t) — u(0,t) remains a 
spectrally negative Levy process for all t > 0. This closure property was first 
noted by Carraro and Duchon in connection with their notion of statistical 
solutions to Burgers equation [11]. That these statistical solutions agree with 
the Cole-Hopf solution for spectrally negative data was shown by Bertoin [5, 
Thm. 2]. The closure property fails if uq has positive jumps — These posi- 
tive jumps open into rarefaction waves for t > 0, and this is incompatible 
with the rigidity of sample paths of Levy processes. An interesting formal 
analysis of closure properties of Burgers equation is presented in [13]. 

Henceforth, we write v(x, t) = u(x, t) — u(0, t) for brevity. The Levy- 
Khintchine representation now implies that the law of the Levy process x i— > 
v(x, t) is completely described by a corresponding "Levy triplet" (bt, o~1, At). 
The mean drift bt = E(u(l, t)) satisfies bt = &o = for every t > 0. Moreover, 
for every t > 0, v(-,t) is of bounded variation, thus the variance <Tj = 0. 
Consequently, the law of v(-,t) is completely determined by only the jump 
measure At which contains the shock statistics. 

It is a striking fact, implicit in [5], that the evolution of A t is described 
by Smoluchowski's coagulation equation with additive kernel, an equation 
that arises in entirely different areas such as the analysis of algorithms [14] , 
the kinetics of polymerization [37], and cloud formation from droplets [23] 
(see [2] for a review). What this means is that mean- field theory is exact 
for Burgers equation with initial data of the form (3), i.e., random one- 
sided data with stationary and independent increments. We give a precise 
statement to this effect below in Theorem 2. Several connections between 
stochastic models of coalescence and Burgers turbulence are reviewed in [8] . 




0, x < 0, 

a spectrally negative Levy process, x > 0, 



(3) 




q > 0, (4) 
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Here, we use the closure property as a basis for a rigorous study of 
universality classes for dynamic scaling in Burgers turbulence. Our motiva- 
tion is the following. A central theme in studies of homogeneous isotropic 
turbulence in incompressible fluids is the universality of the Kolmogorov 
spectrum [30]. A possible rigorous formulation of such universality involves 

(a) the construction of stochastic processes that mimic a 'typical turbu- 
lent flow', and (b) a characterization of the domains of attraction of these 
processes. For Burgers turbulence, step (a) consists of constructing exact 
solutions for special initial data, say white noise or Brownian motion. In 
this article, we carry out step (b) for initial data that satisfy (3). 

Domains of attraction are studied in the classical limit theorems in prob- 
ability (e.g., the central limit theorem), and their process versions (e.g., 
Donsker's invariance principle). For Smoluchowski's coagulation equation 
with additive kernel, we characterized all possible domains of attraction 
in [32], a result akin to the classical limit theorems. In this article we 
deal with a process version. In all that follows, we consider the processes 
x i ^ v(x, t) as elements of the space D of right continuous paths M + — ► R 
with left limits (cadlag paths) equipped with the Skorokhod topology [28, 
Ch. VI]. The shock statistics determine completely the law of this process 
(a probability measure on D). Approach to limiting forms will be phrased 
in terms of weak convergence of probability measures on D. 

Among the initial data we consider, the stable processes are of particular 
importance because of their self-similarity. Let X a ,a £ (1,2] denote the 
stable process with Laplace exponent q a (a = 2 corresponds to Brownian 
motion). The corresponding jump measure A(ds) = s- 1 - 01 ds/T(-a) for a < 
2. There is a one-to-one correspondence between (a) these stable processes, 

(b) statistically self-similar solutions in Burgers turbulence, and (c) self- 
similar solutions to Smoluchowski's coagulation equation. Precisely, this 
works as follows. Let a € (1,2], and let T a denote the first-passage process 
for x i ^ X" + x, i.e., 

Ta=mf{y\X° + y>x}. (5) 

Then for the solution to (1) with uq(x) = X% for x > 0, v(x, t) is statistically 
self-similar, with 

V [x, t) 4 tW-^w := t 1 ^" 1 (xt- 1 '? - , t,x>0, (6) 

where (3 = (a — l)/a. Here = means both processes define the same measure 
on D. The process T a is a pure jump Levy process with Levy measure 
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f a {s) ds, where f a is the number density profile of a self-similar solution to 
Smoluchowski's coagulation equation [32, Sec 6]: 

fa(s) = -J2 [ ~ ' " T{\ + k-kP)suwkp, a€(l,2]. (7) 

k=l 

These solutions are related to classical distributions in probability theory by 
rescaling. If p(s; a, 2 — a) denotes the density of a maximally-skewed Levy 
stable law [18, XVII.7] we have [7, 32] 

f a (s) = S ^ 2 p(s ;a,2-a). (8) 

By the Levy-Ito decomposition [4, Thm 1.1] and (6) we may conclude that 
the magnitudes of shocks in u(-,t) form a Poisson point process valued in 
(0, oo) whose characteristic measure is 

AUds) = t^f a (st 1 -^) ds. (9) 

For 1 < a < 2 the self-similar solutions have algebraic tails, with f a (s) ~ 
s _1-a /r(— a) as s — > oo. The case a = 2 is particularly important since it 
corresponds to Brownian initial data. Here we obtain a solution found by 
Golovin in a model for cloud formation from droplets [23], 

/ 2 ( 5 ) = (4 7 r)- 1 /2 s -3/2 e -s/4_ (1Q) 

For the corresponding solution to (1), the law of v(x,t) can be recovered 
from the law of T%, the first-passage time for Brownian motion with unit 
drift, which is explicitly given as follows (see section 2.6): 

nT*e(y,y + dy)) = ^^exp(-^p¥-) dy. (11) 



Ay 

Considering now arbitrary solutions to (1) with initial data (3), we clas- 
sify solutions that approach self-similar form as t — > oo as follows. We 
establish necessary and sufficient conditions for convergence of the laws of 
the rescaled processes 

X^VP = ^v{\{t)x,t) (12) 

in the sense of weak convergence of measures on D. (Since the shocks co- 
alesce, a rescaling \{t) — > oo is needed to obtain a non-trivial limit. The 
amplitude scaling is natural, see Section 3.) Convergence to a process V* is 

written V® -=> V* as in [28]. Recall that a positive function L is said to be 
slowly varying at oo if lim^oo L(tx)/L(t) = 1 for all x > 0. 
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Theorem 1. Let uq be a spectrally negative Levy process with zero mean 
E(ito(x)), variance ctq > ; and downward jump measure satisfying / °°(s A 
s 2 )A (ds) < oo. 

1. Suppose there is a rescaling X(t) — > oo as t — ► oo and a Levy process V* 
with zero mean E(V]*) such that the random variables converge to 
V* in law. Then there exists a G (1, 2] and a function L slowly varying 
at infinity such that 

al + I r 2 A (dr) ~ s 2 - a L{s) as s -> oo. (13) 
J o 

2. Conversely, assume that there exists a £ (1,2] and a function L slowly 
varying at infinity such that (13) holds. Then there is a strictly in- 
creasing rescaling A(t) — ► oo such that V® — > V a . 

Remark 1. Since and are Levy processes, we have V(*) — > V* if 
and only if we have convergence in law of the random variables V$ for some 
fixed xq £ (0, oo) (see (58)-(59) in section 3 below). We take xq = 1 without 
loss of generality. Part (2) implies in particular that the only possible limits 
are statistically self-similar. 

Remark 2. We say a solution has finite energy if for any finite interval 
I C M.+ we have E (J^ \v(x, t)\ 2 dx) < oo. The jump measure A t for the 
solution is related to the energy by (see Section 4) 

E {^J v(x, t) 2 dx^j = (^J s 2 A t (ds)j j x dx. 

The integral in (13) is thus a measure of the energy in an interval. If it 
is initially finite, it is conserved for t > 0, and it remains infinite if it is 
initially infinite. The only self-similar solution with finite energy corresponds 
to a = 2, and Theorem 1 implies it attracts all solutions with initially 
finite energy. In this sense, one may say that the finite energy solution is 
universal. However, Theorem 1 also indicates the delicate dependence of 
the domains of attraction on the tail behavior of Ao- Heavy-tailed solutions 
seem to us no less interesting than those with finite energy. Finer results 
on asymptotics, and a compactness theorem for subsequential limits that 
builds on Bertoin's Levy-Khintchine classification for eternal solutions to 
Smoluchowski's equation [7], will be developed elsewhere. 
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Remark 3. The case of zero mean, bo = 0, is the most interesting. If 60 > 
or bo < we can reduce to this case by a change of variables (see below). 
If 60 < 0, the solution is defined only for < t < —b^ . Theorem 1 then 
characterizes the approach to self-similarity at the blow-up time. If bo > 
then the behavior of the solution as t — > 00 is determined by the zero-mean 
solution with the same a\ and Ao at the finite time 6q 1 . 

Remark 4. The Cole-Hopf solution is geometric and Theorem 1 may be 
a viewed as a limit theorem for statistics of minima. The utility of regular 
variation in such problems is widely known [33]. If the initial data is white 
noise, the Cole-Hopf solution is a study of the parabolic hull of Brownian 
motion. Groeneboom's work on this problem [24] is the basis for several 
results on Burgers turbulence (in particular [3, 19, 22]). We have been unable 
to find a similar reference to the problem we consider in the probability 
literature ([9] seems the closest). 

Remark 5. There is a growing literature on intermittence, and the asymp- 
totic self-similarity of Burgers turbulence, see for example [21, 26]. Numeri- 
cal simulations and heuristic arguments suggest that this is a subtle problem 
with several distinct regimes. It is hard to obtain rigorous results for gen- 
eral initial data. Theorem 1 tells us that the approach to self-similarity is 
at least as complex as in the classical limit theorems of probability. 

The rest of this article is organized as follows. We explain the mapping 
from Burgers equation to Smoluchowski's coagulation equation in Section 2. 
This is followed by the proof of Theorem 1 in Section 3. Finally, in Sec- 
tion 4 we compute a number of statistics of physical interest: energy and 
dissipation in solutions, the Fourier-Laplace spectrum, and the multifractal 
spectrum. 

2 Mean field theory for Burgers equation 

In this section we explain the connection between Burgers equation with 
spectrally negative Levy process data and Smoluchowski's coagulation equa- 
tion. The main results are due to Bertoin [5] and Carraro and Duchon [12]. 
We follow Bertoin's approach, and explain results implicit in [5] and [7]. We 
think it worthwhile to make this connection widely known in full generality, 
since the results are of interest to many non-probabilists. Exact solutions 
of this simplicity are also useful as benchmark problems for numerical cal- 
culations. 
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2.1 Shock coalescence and Smoluchowski's coagulation equa- 
tion 



Smoluchowski's coagulation equation is a widely used mean-field model of 
cluster growth (see [2, 16] for introductions). We begin with a heuristic 
derivation of the coagulation equation as a mean-field model of shock coa- 
lescence. First consider the evolution of a single shock of size s > 0. Let 
uo(x) = — s1 x >q. Then the solution is 

u(x,t) = -sl x > Xl (t), xi(t) = --t. (14) 

Shock coalescence is nicely seen as follows. Let uq(x) = — Yl k =i s k(0)l x > Xk (o) 
where s k (0) > for k = 1, . . . , N and X\(0) < . . . xn(0). The solution may 
be constructed using the method of characteristics and the standard jump 
condition 

■ l r - +^ 
X = -[U + IT) 

across a shock at x = x(t), where u~ and u + denote respectively the left 
and right limits of u(-,t) at x. At any time t > 0, there are N(i) < N(0) 
shocks at locations x\(t) < Xk(t) < x N ^(t) and 



N k-l S k (t ) 

u(x,t) = - ^2s k (t)l x > Xk( t), x k (t+) = - ^2sj(t+) - k 2 ■ ( 15 ) 

k=l j=l 

The shock sizes s k (t) are constant between collisions, and add upon collision — 
when shocks k and k+1 collide, we set s k (t + ) = s k (t-)+s k+ i(t-) and relabel. 
This yields an appealing sticky particle or ballistic aggregation scenario. We 
say a system of particles with position, mass and velocity (x k (t) , m k (t) , v k (t)) 
undergoes ballistic aggregation if (a) the particles move with constant mass 
and velocity between collisions, and (b) at collisions, the colliding particles 
stick to form a single particle, conserving mass and momentum in the pro- 
cess. We map this shock coalescence problem to a sticky particle system by 
setting m k = s k and v k = x k . Suppose particles k and k + 1 meet at time t. 
Then, with unprimed variables denoting values before collision and primed 
variables denoting values after, since v k+ \ = v k — (m k + m k+ i)/2 we have 

m k v k + m k+ iv k+ i = (m k + m k+1 ) (v k - mk + 1 ^j = m' k v' k . 

The calculations so far involve no randomness. Suppose now that the 
shock sizes Sj are independent and let f(s, t) ds denote the expected number 
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I 




Figure 1: Binary clustering of shocks 

of shocks per unit length with size in [s, s + ds]. We derive a mean-field rate 
equation for / as follows. Let I be an interval of unit length. The number 
density changes because of the flux of shocks entering and leaving / and 
because of shock collisions within /. On average, the velocity difference 
across / is 

poo 

Mi{t)= / sf(s,t)ds, 

JO 

therefore the average influx is Mi(t)f(s,t)ds. Next consider the formation 
of a shock of size si + S2 by a collision of shocks of size si and S2 as shown 
in Figure 1. The relative velocity between these shocks is (s\ + s 2 )/2 (see 
Figure 1). The expected number of neighboring pairs with sizes in [s±, s± + 
dsi], [32, S2 + ds2] respectively is 

f(s 1 ,t)f(s 2 ,t) dsi ds 2 . 

The probability that these neighboring shocks are near enough to collide in 
time dt is \{s\ + S2) dt, thus the number of these shocks that collide in time 
dt is 

f(si, t)f(s 2 , t) Sl 2 ds! ds 2 dt. (16) 

Summing over all collisions that create shocks of size s = s\ + S2, and 
accounting for the loss of shocks of size s (= s± or s 2 ) in collisions with 
other shocks, we obtain the rate equation 

d t f(s,t) = M 1 (t)f + Q(f,f) 

where Q(f,f) denotes the collision operator given by 

2 f'S poo 

Q(f,f)(s,t) = - ^ s/(si,*)/(s-*i,t)<fei- J q (s + s 1 )f(s,t)f(s 1 ,t)ds 1 . 
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We integrate in s to find M\ = Mf, therefore the normalized density //Mi 
satisfies the equation 

sW£) =<3 G4'£)' <17 » 

Up to a change of time scale, this is a fundamental mean-field model of 
coalescence: Smoluchowski's coagulation equation with additive kernel. We 
treat this equation in greater depth below. 

More precisely, it turns out that the random solution u(x,t) has the 
structure described in (15) when the initial data uq consists of a compound 
Poisson process with only downward jumps. The mean drift rate at time t 
is then — Mi(i), and this example shows that the solution blows up at the 
time Mi(O)^ 1 in this case. We show below (see (23)) that one may remove 
the mean drift by a change of scale and slope, yielding 'sawtooth' data 
with a deterministic upward drift that compensates the random downward 
jumps. For such data we obtain a global solution. Thus, there is no essential 
distinction between sawtooth data and the decreasing initial data considered 
above. 



2.2 The Cole-Hopf formula 

The modern notion of an entropy solution stems from the penetrating anal- 
ysis by Hopf of the vanishing viscosity limit to (1). His work was based on a 
change of variables (re)discovered independently by Cole and Hopf [15, 27]. 
This solution is obtained via minimization of the Cole-Hopf function [15, 27] 

H(y,t;x) = i^—^ + j^ U0 (y>) d y'. (18) 

The minimum in y is well-defined for all t > provided U{y) = J'q uo(y')dy' 
is lower semicontinuous and lim x ^± 00 y~ 2 U(y) = 0. This is a mild assump- 
tion and holds for the random data we consider provided that the mean drift 
is zero. We denote the extreme points where H is minimized by 

a_(x, t) = mf{z\H(z, t; x) = min H}, a + (x, t) = sup{z\H(z, t; x) = mini?}. 

y ' y 

(19) 

Notice that any z£R such that x = tuo(z) + z is a critical point of H, and 
represents a Lagrangian point that arrives at x at time t. Of these z, the 
'correct' Lagrangian points are the minimizers of H. If a-(x,t) = a+(x,t), 
this point is unique, and we have 

x — a±(x,t) m , . 

u(x,t) = p ' \ xeR, t>0. (20) 
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There is a shock at (x,t) when a-(x,t) / a + (x,t). In this case, the La- 
grangian interval [a_(x, t), a+(x, t)] is absorbed into the shock and the ve- 
locity of the shock is given by the Rankine-Hugoniot condition (conservation 
of momentum) 

, x u(x+,t)+u(X-,t) 1 fO.+ (x,t) 

u(x, t) = 1 +' ; V 1 = — — -— / no y dy. 21 

2 a+{x,t) -a-{x,t) J a _ {x , t ) 

It will be convenient for us to assume that u is right-continuous and we call 
a(x,t) = a+(x,t) the inverse Lagrangian function. Of course, the speed of 
shocks are still determined by the right-hand side of (21). 

In order to deal with non-zero mean drift in initial data, we will use the 
following interesting invariance of Burgers equation. Assume that uq(x) = 
o(\x\) as \x\ — ► oo, and let u(x,t) be the Cole-Hopf solution with u(x, 0) = 
uo(x), defined for all t > 0. Let c£K and define 

4 c) (x) = u (x) + cx, T c = | I < (22) 

(c) 

Then the Cole-Hopf solution with initial data n is given by 

„(c) (x , i) = _2_„(_^,_^_) + _!_, te|0 ,T c ). (23) 

This is seen as follows. An elementary calculation shows that the Cole-Hopf 
functionals for the different data are related by 



H^(y,t;x) =H(y,- 



X \ cx 2 

+ 



+ ct 7 l + ctj 2(1 + cty 
which implies the inverse Lagrangian functions are related by 

a^(x,t) = a (-?—,—LJ). (24) 
V ' ' \l + ct' l + ctj y ' 

We now substitute in (20) to obtain (23). 

2.3 Solutions with Levy process initial data 

Here we describe how the solution of (1), with initial data of the form (3), 
is determined in terms of Laplace exponents, essentially following Bertoin's 
treatment in [5]. 
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Suppose x i ^ uq(x) is an arbitrary spectrally negative Levy process 
for x > 0, with Laplace exponent ipo having downward jump measure Ao. 
We first show that we may assume without loss of generality that J °°(s A 
s 2 )A()(ds) < oo. Indeed, if sAo(ds) = oo, then uq{x)/x — > — oo almost 
surely as x — > oo. (This follows from the fact that for the compound Poisson 
process X x with jump measure Ao(c2s)1i s i>i, one has X x /x — > oo as x — ► oo 
by the law of large numbers.) In this case the Cole-Hopf function H(y, t; x) 
has no minimum for any t > 0, and equation (1) has no finite entropy 
solution for any positive time. Hence, we may suppose that sAo(ds) < 
oo. 

Next, we show that one may assume the mean drift bo = E(uo(l)) is zero. 
If bo is nonzero, we have lim^oo uq(x)/x = bo a.s. by the strong law of large 
numbers. If &o < 0, then by comparison to compression-wave solutions with 
initial data A + 6max(x,0), we find using the maximum principle that a.s. 
the solution blows up exactly at time — 6q . If bo > there is a global 
solution. In either case, we may use the transformation (23) with c = bo to 
reduce to the case bo = 0, replacing uo(x) by uo(x)—boxl x> o- More precisely, 
we apply (23) for x > noting that a(0, t) > 0, thus a(x, t) > a(0, t) > for 
x > 0, so that (24) holds for x > 0. We have: 

Lemma 1. // is a spectrally negative Levy process with Levy triplet 
(c, <7q, Ao), the Cole-Hopf solution u^(x,t) is determined via (23) for x > 
and t G [0, T c ), in terms of a solution u(x,t) having zero mean drift and 
defined for all t > 0. 

With these reductions, we may restrict ourselves to Laplace exponents 
ipo °f t ne form 

Mq) = J} f+j o {e-i s -l + q S )Ao(ds), q>0. (25) 

We will always assume that a and u are right continuous (i.e., a(x,t) = 
a(x + ,t), compare with (21)). This ensures a is an element of the Skorokhod 
space D, so that we may use the standard Skorokhod topology to study 
limiting behavior. For brevity, we write 

v(x, t) = u(x, t) — u(0, t), l(x, t) = a(x, t) — a(0, t), 

and rewrite (20) as 

v ( Xj t ) = x - l ^ f \ x > o, t>0. (26) 
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Bertoin has shown that for all t > 0, x i— > l[x, t) is an increasing Levy process 
(a subordinator) with the same law as the first passage process for tuo(x)+x. 
We denote the Laplace exponents of / and v by and ip respectively: 

E ^ e -ql(x,t)^ =e -x*(5,t) ) ]g ^ e5 w(x,i)^ = e ^(9,t) j X) t > 0. (27) 

We combine (26) and (27) to obtain 

lKg,t) = f-*(f,f). (28) 
Since / is a subordinator, it has the simpler Levy-Khintchine representation 

/•OO 

<Z>(q,t)=d t q + (l-e-^Mds), q > 0, (29) 
■/ 

where d< > supplies the deterministic part of the drift, and jit is the Levy 
measure of l(-,t), which now must satisfy Jq°°(1 A s)/j,t(ds) < oo [4]. We see 
from (26) that v(-,t) is a Levy process with no Gaussian component, and 
thus has a Levy-Khintchine representation 

ij(q,t)=b t q+ (e- qs - 1 + qs) A t (ds), q>0,t>0, (30) 
J 

related to (29) by 

/•oo ] — d 

A t (ds) = m(tds), h + / sAt(ds) = — — (31) 

Jo t 

Due to the result that I and the first passage process of tuo(x) + x have 
the same law, a simple functional relation holds between i/jq and <&(q,t) [5, 
Thm. 2]: 

MtHg, *)) + *) = ?, ? > 0, t > 0. (32) 

The evolution takes a remarkably simple form when we combine equations 
(28) and (32) to obtain 

i>(q,t)=Mq-tHq,t)), q>0,t>0. (33) 

But then ip(q,t) solves the inviscid Burgers equation (in q and t\) 

d t ip + ipd q ip = o, v(o, q) = Mq)- (34) 

The solution to (34) may be constructed by the method of characteristics 
and takes the form (33). The remarkable fact that the Laplace exponent 
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is also a solution to Burgers equation was first observed by Carraro and 
Duchon [12, Thm. 2]. 

Since ipo is analytic and strictly convex, the solution (33) is analytic for 
all time and unique, and the condition d q ip(0,t) = is preserved for all 
t > 0. By (28)-(31), we have 

poo 

b t = 0, d t = l-t sA t (ds), t > 0. (35) 
J o 

Let 

M = lim Vo(<?), (36) 

q^oo 

We find $(q,t) — ► oo as q — > oo from (32), and differentiate to obtain 

d t = lim d Q <S>(q,t) = lim = - — t > 0, (37) 

with the understanding that dt = if Mq = +oo. Then 



M(t) := lim d q ^{q, t) = [°° sA t (ds) = M ° , t > 0, (38) 

q->oc J Q 1 + tM 



with the understanding that M(t) = t 1 when M = oo. Note M' = — M 2 . 
Below we will characterize the evolution of A t differently. 



2.4 BV regularity 

It is clear from the Cole-Hopf formula that u is locally of bounded variation 
for every t > 0. We derive a decay estimate that quantifies this. The sample 
paths of uq have unbounded variation if and only if [4, p. 15] 

<rj>0 or / sA (ds) = oo. (39) 
Jo 

Heuristically, this corresponds to the presence of many small jumps ('dust'). 
This is reflected in the Laplace exponent as Mq = lim^oo vp' (q) = +oo in 
this case. On the other hand, Mo is finite if and only if no is BV, in which 
case (To = and Mo = J °° sAo(ds) < oo. 

The analytic formula (29) has the following probabilistic meaning. If 
we take a Poisson point process x *— > m* (masses of clusters) with jump 
measure we have the representation [4, p. 16] 

l(x,t) = d t x + ^2 mt y ( 40 ) 

0<y<x 
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The velocity field, and a point process of shock strengths s* = t 1 m t y are 
determined from (20), (37) and (40) by 

v(x,t) = M(t)x-j m* =M(t)x- J^s*, (41) 

0<y<x 0<y 

For every t > 0, v(x, t) is the difference of two increasing functions: a 
linear drift and a pure jump process. Thus, it is of bounded variation, and 
by (35) and (40) we have 

eQT \d y v(y,t)\dy^ =2M(t)x, x,t > 0, (42) 

because E (j2o<y<x m y) = xM(t) = x(l - dt). 

2.5 Relation to Smoluchowski's coagulation equation 

We consider a positive measure v T {ds) interpreted as the number of clusters 
of mass or size s per unit volume at time r. Clusters of mass r and s coalesce 
by binary collisions at a rate governed by a symmetric kernel K(r,s). A 
weak formulation of Smoluchowski's coagulation equation can be based on 
a general moment identity for suitable test functions £ (see [32]): 

C(s)u T (ds) = 

-| poc roc 

~ / / (C(r + s) ~ C(r) ~ COO) K(r, s) v T (dr) v T {ds). (43) 
1 Jo Jo 

We consider only the additive kernel K(r, s) = r + s. It is classical that 
(43) can then be solved by the Laplace transform [16]. We denote the initial 
time by To (to be chosen below). The minimal (and natural) hypothesis on 
initial data v T0 is that the mass J Q sv T0 (ds) is finite. We scale the initial 
data such that J °° sv T0 = 1. The Laplace exponent 



<p(q,r)= / (l-e-«> T (da) (44) 
Jo 

then satisfies 

d T <p - Lpdgif = -Lp, t > Tq. (45) 

We showed in [32] that (45) may be used to define unique, global, mass- 
preserving solutions to (43). In particular, a map r i— > v T from [tq, oo) to 
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the space of positive Radon measures on (0, oo), such that J °° sv T {ds) = 1 
for all r > To, is a solution of Smoluchowski's equation with K(r, s) = r + s 
(in an appropriate weak sense detailed in [32]) if and only if ip satisfies (45). 

We now connect solutions of the inviscid Burgers equation (1) with Levy 
process initial data to solutions of Smoluchowski's equation through a change 
of scale. Let uq satisfy (3) and assume as in subsection 2.3 that the corre- 
sponding downward jump measure Ao satisfies f^°{s A s 2 )Aq(cIs) < oo and 
the mean drift is zero. Let At be the jump measure of the Cole-Hopf solu- 
tion. With M and M(t) as in (36) and (38), let r = -logM if u is of 
bounded variation, and To = — oo otherwise, and set 

r = -logM(t), u T {ds) = A t (M(t)ds). (46) 

From (38) it follows / °° sv T (ds) = 1, and by (44) and (30) we find 

V (q,r)=q-^(qe T ,t). (47) 

We see that ip solves (34) if and only if ip solves (45). Therefore, the rescaled 
Levy measure of v(-,t) evolves according to Smoluchowski's equation. Con- 
versely, given any solution of Smoluchowski's equation with initial data uq 
at a finite To, we can construct a corresponding solution of (1) by choosing 
uq to be a spectrally negative Levy process with jump measure A to via (46) 

Initial data uq with unbounded variation are of particular interest. Here 
we have eternal solutions v T to (43) defined for all r£l. We see that eternal 
solutions are in one-to-one correspondence with initial data uq of unbounded 
variation via (47). A finer correspondence mapping the clustering of shocks 
to the additive coalescent is found in [6]. 

To summarize, we have the following correspondence. 

Theorem 2. Assume uq is a spectrally negative Levy process with Levy 
triplet (0,<7q,Ao), with the same assumptions as in Theorem 1. Then for 
all t > 0, v(-,t) and l(-,t) are Levy processes with triplet (0,0, At), whose 
jump measures A t determine a solution v T {ds) to Smoluchowski's coagulation 
equation with rate kernel K(r, s) = r + s as described in (46). 

2.6 Self-similar solutions 

Bertoin's characterization of eternal solutions is the analogue of the Levy- 
Khintchine characterization of infinitely divisible distributions [4, 18]. Among 
the latter, the stable distributions are of particular interest, and their ana- 
logues for Smoluchowski's equations are obtained by choosing the Laplace 
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exponent ipo(q) = q a , a G (1,2]. For a G (1,2) the corresponding Levy 
measures are 

a -(i+«) 

A(ds) = — da. 

r(-a) 

The Laplace exponent q 2 corresponds to an atom at the origin. We thereby 
obtain for a G (1,2] a family of self-similar solutions to Smoluchowski's 
equation with Laplace exponent of the form (p(r,q) = e~ l3T (p a (qe l3T ) where 
(p a solves 

<p a (q) a + <p a (q) = q, q>0- (48) 
The self-similar solutions to Smoluchowski's coagulation equation are 

u T (d S ) = e- 2T ^Ue- T ^s)ds, /3 = — , a€(l,2], (49) 

a 

where / Q has been defined in (7). An analytic proof that these are the 
only self-similar solutions to Smoluchowski's equation may be found in [32]. 
Each of these solutions corresponds to a self-similar process. Precisely, let 
X a denote the stable process with Laplace exponent q a , and T a and V a 
denote the processes 

T? = w£{y>0:XZ + y>x}, V x a = x-T«. (50) 

We have Mo = +00 and M(t) = t" 1 = e~ T in this case, and the Laplace 
exponent of the process l(-,t) is of the self-similar form 

^{q,t) = ^{q,T)=r 1 '^ a (qt 1 ^), t > 0. (51) 

The solution processes have the scaling property 

l(x,t) ± tWT«_ 1/g , v(x,t) k tW-iV^. (52) 

The corresponding Levy measures are obtained from (31), (46) and (49): 

tf(ds) = t-Wf a {t-Ws) ds, AUds) = t 1 - 2 ^^ (f-Ws) ds. (53) 



In the important case a = 2, we have 1//? = 2 and f2(q) = —\ + y \ + q, 
and by Laplace inversion [1, Ch.29] we obtain the explicit expression in (11) 
for the distribution of T 2 . 
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3 The convergence theorem 

In [32] we proved the following theorem characterizing solutions that ap- 
proach self-similar form in Smoluchowski's coagulation equation with addi- 
tive kernel. To every solution v T of (43) with J °° su T {ds) = 1 we associate 
the probability distribution function 

F(s,t) = / rv T (dr). (54) 

J(0,s] 

To a self-similar solution f a , a G (1, 2] with (3 = (a — l)/a we associate 

F a (s) = [ rf a {r) dr = f j ^^ (1 + k - (55) 

A probability distribution function F* is called nontrivial if F*(s) < 1 for 
some s > 0; this means the distribution is proper (lim^oo F*(s) = 1) and 
not concentrated at 0. 

Theorem 3. Suppose t± G R and v T , r G [ri,oo) ; is a solution to Smolu- 
chowski's coagulation equation with additive kernel such that f °° sv T1 (ds) = 
1. 

L Suppose there is a rescaling function A(r) — > oo as r — > oo and a 
nontrivial probability distribution function F* such that 

lim F(X(t)s,t) =F*(s) (56) 

T— »00 

af aZZ points of continuity of F* . Then there exists a E (1,2] and a 
function L slowly varying at infinity such that 

r 2 u T1 (dr) ~ s 2 ~ a L(s) as s — ► oo. (57) 

2. Conversely, assume that there exists a G (1,2] and a function L slowly 
varying at infinity such that (57) holds. Then there is a strictly in- 
creasing rescaling A(r) — ► oo such that 

lim F(X(t)s,t) = F a (s), < s < oo, 

T— >00 

where F a is a distribution function for a self similar solution as in 
(55). 



L 
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We now prove Theorem 1. Let uo be a spectrally negative Levy process 
with zero mean drift and J °°(s A s 2 )Ao(ds) < oo. To the solution increment 
v(x,t) = u(x, t) — u(0, t) with downward jump measure At, associate a solu- 
tion v T of Smoluchowski's coagulation equation (43) as in Theorem 2 with 
Laplace exponent <£>(<?, t) given by (47). Let t\ = To = — logM if M < oo, 
and let n = if Mq = +oo and To = — oo. 

We deduce Theorem 1 from Theorem 3 by establishing two equivalences: 

(a) There is a rescaling A(i) — > oo as t — > oo and a Levy process V* with 

zero mean drift E(V*) such that — > V* if and only if there is a 
rescaling A(r) — > oo as t — > oo and a nontrivial probability distribution 
function F* such that (56) holds. 

(b) J °° s 2 z/ Tl (ds) < oo if and only if J °° s 2 Ao(ds) < oo. Moreover, f£ r 2 u n (dr) ~ 
s a L(s) as s — > oo if and only if J Q S r 2 Ao((ir) ~ s 2 ~ a L(s) as s — > oo. 

Proof of (a). We prove claim (a) by showing each part equivalent to a 
corresponding convergence statement for rescaled Laplace exponents. First, 
convergence in law in B for processes with independent increments can be 
reduced to the convergence of characteristic exponents [28, Cor. VII. 4. 43, 
p. 440]. In particular, suppose A(i) — ► oo as t — ► oo. Then we have 

v® 4. v*, with E(y;) = o, (58) 

if and only if E,(e tkV * ) — ► E,(e tkVx ) for all fc£R, uniformly for x in compact 
sets, and E(VJ*) = 0. But since we are working with Levy processes, the 
Levy-Khintchine formula shows the dependence on x is trivial, and thus (58) 
is equivalent to 

E(e ifcV i (t) ) -► E(e ifey i*) for all keR, and E(F X *) = 0. (59) 

But pointwise convergence of characteristic functions is equivalent to con- 
vergence in distribution of the random variables vj^ [18, XV. 3. 2], and since 
Vl t] = 1 - T-P < 1, (59) is equivalent to convergence of the Laplace trans- 
forms [18, XIII. 1.2]: 

E(e« y i (t) ) -» E(e qV i) for all q > 0, and E(V X *) = 0. (60) 
Taking logarithms and using (12) and (27), (60) is equivalent to 

A^ (qt/\,t) ^>*(g) for all q > 0, and d g V*(0) = 0, (61) 
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where R(e qV * ) = e^*^. This expresses the convergence of V® in terms of 
convergence of rescaled Laplace exponents. 

Now suppose A(t) — > oo as r — > oo. Using [18, XIII. 1.2] again, the 
(proper) convergence in (56) is equivalent to pointwise convergence of Laplace 
transforms: 

rj(q,T) -► rf(q) for all g > 0, with rj*(0) = 1. (62) 

where r/(g,T) := J °° e - gs F(A(r) ds, r), r/*(g) := J °° e -« s F*(ds). By (54) 
and (44), we have 

rj(q,T) = (d q tp)(q/\,T), [ r](r,T)dr = \tp(q/\,T). (63) 

JO 

We claim that (62) is equivalent to the statement that (with <p*(q) = 
foV*{r)dr) 

X(p(q/X,r) (p*(q) for all q > 0, and 0,p*(O) = 1. (64) 

Clearly, since ??(-, r) is completely monotone and bounded, (62) implies (64). 
In the other direction, assume (64). For any sequence Tj — > oo there is 
a subsequence along which rj(q,Tj) converges for all (rational, hence real) 
q > 0, to some limit whose integral must be p* . Thus (62) follows. 

We now finish the proof of claim (a) by observing that due to (47), we 
have 

\<p(q/\,T)=q-W(qe T /\,t). (65) 
Hence the convergence in (61) is equivalent to that in (64) provided we have 

\(t)/t = \(r)/e\ (66) 

or \{t) = tM(t)X(r), since tM(t) -> 1 as t -> oo. (Note tM(t) = 1 if 
M = oo.) 

Proof of (b). It is only the case M$ = oo that requires some work. Indeed, 
if Mq < oo we see from (46) that v Tl (ds) = Ao(Mo ds). In what follows, we 
suppose that Mq = oo. We then have an eternal solution to Smoluchowski's 
equation, and t = e T . We shall compare the tails of uq (r = 0) with that of 
Ao (/ = 0). 

Claim (b) is a purely analytic fact that follows from Karamata's Taube- 
rian theorem [18]. We first reformulate it as a statement about Laplace 
transforms. Let tfo(q) = <p(q, 0), ipoio) = ^(<2S0). For every a £ (1,2] we 
have 

f r'Mdr) ~ s 2 ~ a L(s) <=> 1 - ^(g) ~ (-) (67) 
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as s —> oo and q — > respectively (see [32, eq. 7.4]). By the same argument, 
f r 2 A (dr) ~ s 2 -°L( S ) ^ ~ (f^L (-) F ( 3 ~ ") , (68) 

with the following caveat when a = 2. If f °° r 2 A(dr) = oo then (68) holds. 
On the other hand, if / °° r 2 Ao(dr) < oo then we must modify the second 
condition in (68) to 

tfo(q) ~ (o 2 + J r 2 A(dr)^jq, q - 0. 
We set t = 1 in (47) and differentiate (33) with respect to q to obtain 

MM*)) = 1 ~^{ q) = -7TT " 1- ( 69 ) 
fPoW) VW) 

The functions ^>q, (fo , and are strictly increasing. Since v?q(0) = 1 we 
also have <po(q) = <?(1 + o(l)) as q — ► 0. A sandwich argument as in [18] may 
now be used to deduce claim (b). First suppose that (67) holds. Fix b, e > 0. 
Then for q sufficiently small we use monotonicity and (69) to obtain 

l-^{bq{l-e)) <f/ (q(l+e)) < < 1 ~ + *)) <So(q(l ~ *)) 



1 - ^(9(1 + e)) <p> (bq(l - e)) %{q) 1 - ^(g(l - e)) ^(6,(1 + e)) ' 
Letting first q and then £ — ► 0, we obtain 

]im jmi =b a-l 
g^O i// Q (q) 

Thus, ip' is regularly varying with exponent a — 1. Similarly, if we assume 
that (68) holds, we sandwich 

i-pW yd(g) . W(l+£)g) 

^((1 + 1 - tp' (q) <p' (bq) ' ^ ((1 - ' 

to deduce that 1 — is regularly varying with exponent a — 1. Finally, 
since </?q(0) = 1 it follows from (69) that lim g _ > oV'o (<?)/(! — V'dfa)) = 1- This 
finishes the proof of Theorem 1. 



4 Energy, dissipation and spectra 

In this section, we compute several statistics of physical interest for the 
solution increments: mean energy and dissipation, the law of the Fourier- 
Laplace transform, and the multifractal spectrum. While the computations 
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are routine, some interesting features emerge, namely (i) conservation of 
energy despite dissipation at shocks, (ii) a simple evolution rule for the 
Fourier-Laplace spectrum, and (hi) a multifractal spectrum in sharp variance 
with that of fully developed turbulence. Simple proofs of fine regularity 
properties (e.g., Hausdorff dimension of the set of Lagrangian regular points) 
may be found in [5]. 



4.1 Energy and dissipation 

The energy in any finite interval I C M+is computed using the Levy- 
Khintchine formula (30) and Fubini's theorem as follows. 



dx 

9=0/ 



E N v(x, t) 2 dx\ = Je (v(x, t) 2 ) dx = jf (d 2 E 
= / d 2 e x '^ dx= I (x 2 d q iP(0,t) 2 +xd 2 M,t)) dx 

Jl 9=0 J J 

= b 2 J^x 2 dx+^ y 2 A t (dy)^j j^xdx. (70) 

Let us restrict attention to solutions of mean zero, that is b t = 0. Then we 
have conservation of energy in the sense that 

E (^J v(x,t) 2 dx^j = E (^J v(x,0) 2 dx^j , t > 0. (71) 

Indeed, by (70), we see that (71) is equivalent to 

POO 

d 2 M t) = dfyo = al+ s 2 A (ds) =: M 2 , (72) 

Jo 

with the understanding that d 2 ip(0,t) = oo if / °° s 2 Ao(ds) is divergent. It 
is only necessary to differentiate (33) to obtain 

OMq, t) = MS-Z^l d 2 Mq, t) = =Mq - Hp), 

and then take the limit q — > to obtain (71). 

The dissipation at a shock with left and right limits u± is obtained as 
follows. The decay of the L 2 norm for solutions to Burgers equations with 
viscosity e, ut + uu x = eu xx , is given by 

^- f u 2 dx = 2e f u 2 dx. 
dt Jr 
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The right hand side may be evaluated exactly for traveling waves (viscous 
shocks) of the form u(x, t) = u £ {x — ct). It is easily seen that for any e > 
a traveling wave profile connecting the states u_ > u + at =foo is of the form 
u £ (x — ct) = w{(x — ct)/e) where w satisfies the ordinary differential equation 

1 , 2 2 \ d"W U_ + U + 

-c{w-u-) + -(w -u_) =—, c= . 

We therefore have 

2e f u\ dx = 2 f (w'f d£ = 2 f + -c{w - u_) + \ (w 2 - u 2 _) 

Jr ' Jr Ju- l ^ 



dw 



3 



= 2(n_ - U+ ) 3 «;(1 - «;) dw = (U " Q U+) ■ 
Jo 3 

The right hand side is independent of e and captures the dissipation of the 
entropy solution in the limit e — ► 0. The dissipation at shocks in any finite 
interval I C M+ may now be computed by summing over all shocks in / 
using (40) and (41): 

(73) 

where |/| is the length of /. 

Conservation of energy in the sense described in (71) is rather surprising 
in view of the dissipation at shocks. In particular, there are solutions with 
finite energy (J °° s 2 A t (ds) < oo), but infinite dissipation (J °° s 3 A t (ds) = 
oo). However, there is no contradiction, since (71) refers to the expected 
value of the energy in any finite interval /, and the energy dissipated in 
shocks is compensated by energy input from the endpoints of /. 



4.2 The Fourier-Laplace spectrum 

We show that the law of the Fourier transform v(k,t) of paths x ^ v(x,t), 
is determined by a Levy process with jump measure s~ 1 A t (s)ds, where 
A-t(s) = J* s °° At(ds). Here At(ds) denotes the jump measure of v(x,t) (see 
Theorem 4 below). The assertion v ~ k" 1 as k — > oo for white noise ini- 
tial data is common in the Burgers turbulence literature (e.g., see [20, 36]). 
For the present case of Levy process initial data, we show that v(k, t) ~ 
—iM(t)k~ 2 as k — > oo. In addition, we find precise corrections under addi- 
tional assumptions on A t (for example, for self-similar solutions). 



23 



These computations with the laws of the Fourier-Laplace transform should 
be contrasted with the conventional notion of the power spectrum. Despite 
its widespread use for wide-sense stationary processes, the power spectrum 
is of limited utility for the present problem involving stationary increments, 
as we now show. Fix L > and consider the interval [0, L]. Almost every 
sample path v(x,t) is bounded on [0, L] and we may define the truncated 
Fourier transform 

v L (k,t) = f e- ikx v(x,t)dx. (74) 
Jo 

If the energy is finite (M 2 < 00 in (72)), we may compute a truncated power 
spectral density Sx(fc) as follows. We have 

L\v L (k,t)\ 2 = y [ [ e- ik( - x - y ^v(x,t)v(y,t)dxdy. (75) 
L L Jo Jq 

Since v(x,t) is a Levy process with mean zero, the autocorrelation is 

E(v(x,t)v(y,t)) = (xAy)M 2 . (76) 
We take expectations in (75) to find 

SL{k):=^E(\v L {k,t)n = ^(l-^y k^O. (77) 

The power spectrum S(k) = lim^oo Sl(A;) = 2M2/A; 2 is now seen to be well- 
defined, but is unsuitable for distinguishing solutions because all solutions 
with the same energy (possibly infinite) have identical power spectrum. 

A well-defined spectrum that distinguishes solutions may be obtained by 
taking the Fourier-Laplace transform of process paths. For fixed p > we 
define the random variable 



POO 1 /'OO 

Cv(p,t)= e~ px v(x,t)dx = - e- px v(dx,t). 
Jo P Jo 



(78) 



The integrals are well-defined because lim x -, 00 v(x,t)/x = a.s. by the 
strong law of large numbers, and v(x, t) is of bounded variation. If s* denotes 
a point process of shock strengths as in (41) we have 

p£v{p,t) = ^--J2e~ px 4- (79) 

P 0<x 
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We determine the law of pCv(p,t) by computing its Laplace transform via 
the 'infinitesimal' Levy-Khintchine formula E i y e qv( - dx ^ = e ^ q ^ dx . We then 
have 



E (e qpCv ^ = exp °° ^(qe~ px , t) dx 

= exp Q £ ^ =: exp Qv#(<7, *)) . P, 9 > 0, (80) 

after the change of variables </ = qe~ px . We now observe that ip# determines 
a Laplace exponent as follows. Let A t (s) = J s °° A t (ds) denote the tail of the 
Levy measure A t . Since f °°(s A s 2 )A t (ds) < oo we have the bounds 

poo pe 

sA t (s)< rA t (dr), s 2 A t (s) < / r 2 A t (dr) + s 2 A t (e), s€(0,e). 

is JO 



Therefore, 



lim sA t (s) = 0, lim s A t (s) = 0, 

s^oo s^O 



and we may integrate by parts in (30) to obtain 



q> 



poo 

/ {l-e- q ' s )A t {s)ds. (81) 
Jo 

Integrating once more in q' we find 

Jo Q Jo s 

We integrate by parts in (38) to see that 

poo poo 

/ A t (s) ds= sA t (ds) = M(t) < oo. 
Jo Jo 



(83) 



This enables us to write 



f°° A*(s) 
^(q,t) = M(t)q-$ # (q,t), *#(q,t) = (1 - e~ qs )-^ ds. (84) 

Jo s 

Since (83) ensures s~ 1 A t (s)ds satisfies the finiteness conditions for a jump 
measure, ip# is a Laplace exponent for a Levy process with zero mean drift 
that we denote by Z t . Similarly, <3?# is the Laplace exponent for a subordi- 
nator that we denote Y t . We summarize our calculations in the identities 

Z\ = M(t)r-Y?, E (> z ') = e ^#fet), E (e _?1 ?) = e^**^, r,q,t> 0. 

(85) 
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The result is that the Laplace spectrum of the solution increments is 
determined by 



E 



e qpCv(p,t)\ _ jg 




q > 0, p > 



(86) 



which implies that £v(p,t) has the same law as p 1 Z\/ p f° r fixed p > 0. 
Note that Cv(p, t) is not a Levy process in p. In fact, for a fixed realization, 
Cv(p,t) is analytic in p. Nevertheless, its law is determined by the Levy 
process Z l . 

We extend this computation to the Fourier spectrum (p = ik) as follows. 
The calculations leading to (80) hold for complex q with Re(g) > 0, and in 
particular for q = i£, £ £ R. Moreover, j£v(p,t) is a well-defined random 
variable for every p with Re(p) > 0. Thus, we may analytically continue the 
identity E(e q P Cv ^) = exp(p _1 V>#(9, *)) to a11 P with Re (p) > 0, and q = 
i£. As in (4), let t) = —il)#(i£,t) define the characteristic exponent 

corresponding to the Levy process Z l . For e, k > we set p = e + ik, 
v(k — ie, t) = Cv{e + ik, t) and pass to the limit e^Oon both sides of (80) 
to obtain 



Thus, for fixed k > 0, as e [ the random variables ikv(k — ie,t) converge 
in law to the (real) random variable Z*, k . We denote this limit by ikv(k,t). 
As before we do not assert that the processes ikv(k — ie,t) converge in law 
to the process Zy k , simply the convergence of random variables for fixed k. 
We summarize our conclusions as follows. 

Theorem 4. Let Y l be a subordinator with Laplace exponent 3>#(q,t) from 
(84), and let Z l the Levy process defined by (85). Then for every fixed p > 
and k > the random variables pCv(p,t) and ikv(k,t) have the same law 
as Zy p and Zy k , respectively. 

Due to this result and (85), we always have the upper bound ik 2 v(k, t) < 
M(t) a.s. This crude bound may be refined as k — > oo using information 
related to the sample path behavior of subordinators (see [4, Ch. III.4]). 

Corollary 1. For every t > 0, lim^oo ik 2 v(k, t) = M(t) in probability. 



limE 

slO 




(87) 



(88) 
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Proof. This follows from the fact that lim r jo Y* /r = in probability, proved 
as follows. By (85) we have K(e~ qY ^ r ) = e~ r ®*( q l r,t \ and since 1 — e~ s < 
1 A s, by (84) we have that the Laplace exponent 

poo A ( o) p°° A (a\ 

r$ # (q/r,t)=r (1 - e - qs/r )-^- ds < / (rAqs)^^ds^O 
Jo s Jo s 

as r I for each q > 0. Hence Y*/r — > in law. □ 

A similar conclusion holds for the Laplace spectrum as p — > oo. Actually, 
for the subordinator Y r ', the sample paths have the stronger property that 
lim r |oY r < / r — > a.s. [4, III. 4. 8]. Under a mild assumption on the integra- 
bility of the small jumps, we can strengthen convergence in probability to 
almost-sure convergence of the Laplace spectrum. 

Corollary 2. For every t > 0, linip^oo p 2 Cv(p,t) = M(t) in probability. If 
we also assume J \ log s|A t (s) ds < oo, then linip^oop 2 Lv (p, t) = M(t) a.s. 

Proof. For notational convenience, we suppress the dependence on t in the 
proof. Fix e > 0, and let p m = 2 m for positive integers m. We will show 
that \\m m ^ OQ p 2 n Cv{p m ) = M a.s. That is, for every e > 0, we claim 

P (\p 2 m Cv{p m ) - M\ > e infinitely often) = 0. (89) 

This is sufficient to establish linip^oop 2 £v(p) = M a.s. Indeed, since M/p — 
pCv(p) is completely monotone by (79), for p G (p m ,p m +i) we have the 
bounds 

< M(t) - p 2 Cv{p) < V- (M(t) - p 2 m Cv( Pm )) <2(M-p 2 m Cv{ Pm )) , 

Pm 

and therefore 

{M - p 2 Cv{p) > 2e] C{M- P 2 m Cv{p m ) > e}. (90) 
In order to prove (89) we use the elementary estimate 

P (\p 2 m Cv(p) - M\ > e) = P { Pm Y 1/Pm >e)< ^-E (l - exp (-^Y 1/Pm )) 

-M 1 -" (-^* ( t)))^>^ 

We will show that Y^=iPm~^#(Pm/ e ) < 00 • ^he first Borel-Cantelli lemma 
then implies (89). For clarity, we suppose e = 1. This causes no essential 
difference and reveals the main computation. 
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Denote the integrated tail of the Levy measure for Y l by 

n*(-)= r^p-ds'. (9i) 

is s 

We integrate by parts and use Tonelli's theorem to find 

OO .qo oo 

£ PAiPm) = / £ e" PmS n t ( S ) ds. (92) 

1 «^ 1 



m=l J ° m=l 



It is only necessary to check that the integral over s G (0, 1) is finite. Here 
we use the elementary estimate 



°° /'OO -I /» 



10 | logs) + 1 

e y — < 



log 2 ' 
so that 

/ y2e~P™ s Il t (s)ds<— (l + |loga|)nt(a)ds. 
Jo m=1 iog2J 

By the definition of Ht(s) in (91), the last integral is 

/ I log s I / ^—^-drds = [ ^—^-dr [ I log s I ds 
Jo Js r Jo r Jo 

r 1 roc 

< / \log r\A t (r) dr + / A t (r) dr, 
Jo Jo 

which is finite by assumption. □ 

Corrections to the bound ik 2 v(k, t) < M(t) involve the law of the iterated 
logarithm [4, III. 4]. The following corollary holds for initial data that is not 
BV (so Mo = +oo and M(t) = 1/t) with suitably regular small jumps 
('dust'). 

Corollary 3. Assume <ro / and a = 2, or assume ctq = and Ao(s) = 
J* s °° A(dr) is regularly varying at zero with exponent —a where a G (1,2). 
Then for every c > and t > we have 

m f t^ 1 — ik 2 v(k,t) \ log log A; , . . 

— log P [ — ttt~\ — ; ,\ < c Si , k^oo, 93 

B V h(k log log A;) - / 7cTt 1+2 'r' v 7 

where 7 = 1/(q — 1) and h{k) = k/tpo(k)- 
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This corollary is a consequence of [4, Lemma III. 12] and is associated 
with the following lemma of independent interest which shows that the evo- 
lution preserves the regularity of the dust. 

Lemma 2. (a) Assume that uq = and Ao(s) is regularly varying at zero 
with exponent —a, a G (1,2). Then A t (s) is regularly varying at zero with 
exponent — 1/a for every t > 0. 

(b) If <7 / 0, then A t (s) ~ (cr t) _1 y/2/(irs) as s — > 0, for every t > 0. 

Proof. Recall that i _1 — ik 2 v(k,t) agrees in law with kYy k . Combining (28) 
with (82) we find that the Laplace exponent of the subordinator Y l satisfies 

We claim that <&(-,t), and hence $#(-,t), is regularly varying at oo with 
exponent a = 1/a G [5, 1). 

To prove the claim, we integrate by parts in (25) to obtain 



^o(g) 
q 2 



^2 + L°° e ~ qs (/°° Ao(r)dr ) ds ' 



First assume uo = and Ao is regularly varying at zero with exponent 
—a. Then ipo is regularly varying at infinity with exponent a. This fol- 
lows from [18, XIII. 5. 3], or may be proved directly. If ctq 7^ 0, we have 
liniq^oo ipo(q)/q 2 = °o/2- The Laplace exponent $>(q,t) is determined via 
the functional relation (32). Since a > 1, the map i-> <7o($) : = V'o(^) + <£ 
is regularly varying (in <&) at 00 with exponent a. Therefore, the inverse 
function &(q, t) is regularly varying (in q) at 00 with exponent 1/a. 

Now let g#(-,t) be the inverse function to <]?#(•, i). Then by [4, Lemma 
III.4.12] we infer that for every c > 0, 

/ fry 1 \ 
- lo g P I n 1 T i n < c ) ~ (l-d)(d/c) a /( 1_d hoglogfc, k^oo, 



(95) 



h # (klog log k,t) 
where 

By (94) and regular variation we have $#(q,t) ~ <&(q/t,t)/a as g — ► 00, and 
thus by (32) we find that as q — > 00, 

##(<?,*) ~ *5o(%) ~ t^o(atq) ~ t 1+Q a _Q Vo(g)- 
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Substituting c = ct 1+a a~ a into (95) yields Corollary 3. 

Karamata's Tauberian theorem and the monotone density theorem now 
imply that At(s) is regularly varying at zero with exponent —1/a. If ctq 7^ 0, 
we find $(q,t) ~ (aot)- 1 ^ as q — > 00. Assertion (b) of the Lemma then 
follows from the Tauberian theorem. □ 

For the self-similar solutions, ipo(q) = q a with a £ (1,2], Aq(s) = 
S - a /(aT(-a)) for a G (1,2), and we have A t (s) ~ t- 1 (ts)~ 1 / a /r(l - 1/a) 
as s — > 0. 

4.3 The multifractal spectrum 

The notion of a multifractal spectrum was introduced by Frisch and Parisi 
to describe the intermittency of velocity fields in fully developed turbu- 
lence [20]. The multifractal spectrum d(h) measures the dimension of the 
set Sh where the velocity field has singularities of order h. There are differ- 
ent mathematical formulations of multifractality, corresponding to different 
notions of what one means by singularities of order h. Here we follow the 
treatment by Jaffard, which yields d(h) rather easily [29] (the notation has 
been changed slightly for consistency with this article). 

We say a function v : R+ — > R, is C r (xo) for a point xo £ R+ if there is 
a polynomial P XQ of degree at most [r] such that 

\v(x) - P Xo {x)\ < C\x - x | r , 

in a neighborhood of xq. The Holder exponent of v at xq is defined as 

h v (x ) = sup{r \v £ C r (x ) }. 

We define Sh to be the set of points where v is of Holder exponent h. The 
multifractal spectrum d{h) is the Hausdorff dimension of Sh- If Sh is empty, 
the convention is d(h) = —00. As an example, let us compute the multifrac- 
tal spectrum when the initial data is of bounded variation. Then Mq < 00 
in (41) and there is a finite number of shocks s* in a finite interval [0, x] 
with probability 1. Suppose xq is not a shock location for v(-, t). Then (41) 
shows v is analytic near xo and h v (xo) = 00. If xo is a shock location, then 
v G C~ e (xo) for every e > 0, so that h v (xo) = 0. Thus, we have simply 
d(0) = and d(h) = —00 for every h 7^ 0. 

The multifractal spectrum is more interesting for initial data of un- 
bounded variation, that is, when (39) holds. In this case, the jumps in 
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v are dense. Following Jaffard [29], the multifractal spectrum is computed 
as follows. We define 

Cj(t) = C A t (ds), (3 t = max f 0, limsup ] . (96) 

J2-3- 1 \ J log 2 J 

For any t > 0, v(-,t) has no Brownian component. It then follows from [29, 
Thm. 1] that 

d t {h) = {~^ fcepWA], (97) 

Experiments suggest that the multifractal spectrum for fully developed three- 
dimensional turbulence is a concave curve [31, Fig. 2]. This is in clear 
contrast with (97). 

For example, let us compute the multifractal spectrum for the self-similar 
process V a of index a £ (1,2]. Since A°(s) is a scaled copy of Af(s), dt(h) 
is independent of t. We use (7) to obtain the asymptotics as s — > 0: 

A?(cfa) = f a (s) ds ~ ^l s ^T{2 - p) ds, (3 = — . 

7r a 

We then have (3 t = a -1 , i > and 

d(h) = \ h J a > h *M> (98) 

In particular, (98) implies that d{a) = 1, that is v(x,t) is C a (x) for a.e 
x £ R+. For this set a finer characterization of the local variation of v(-,t) 
may be obtained by using the Fristedt-Pruitt law of the iterated logarithm 
(sec [5, Cor. 1]). However, the multifractal spectrum, also describes sets Sh, 
< h < a, that are not covered by the Fristedt-Pruitt law. 
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